TWO-PHASE FLOW PROBLEM COUPLED 
WITH MEAN CURVATURE FLOW 
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Abstract. We prove the existence of generalized solution for incompressible and viscous non- 
Newtonian two-phase fluid flow for spatial dimension 2 and 3. The phase boundary moves along 
with the fluid flow plus its mean curvature while exerting surface tension force to the fluid. An 
approximation scheme combining the Galerkin method and the phase fleld method is adopted. 



1. Introduction 

In this paper we prove existence results for a problem on incompressible viscous two-phase fluid 
flow in the torus n = T'^ = (M/Z)"', d = 2, 3. A freely moving (d— l)-dimensional phase boundary 
r(i) separates the domain 0, into two domains and (t), t > 0. The fluid flow is described 

by means of the velocity field u : $7 x [0, oo) — t- and the pressure 11 : x [0, oc) — t- M. We assume 
the stress tensor of the fluids is of the form r^(u, 11) = r^(e(n)) — 11/ on ft^{t), respectively. Here 
e{u) is the symmetric part of the velocity gradient Vu, i.e. e(n) = {\/u + \7u'^)/2 and / is the dxd 
identity matrix. Let §i{d) be the set oi d x d symmetric matrices. We assume that the functions 
: S{d) — )• §{d) is locally Lipschitz and satisfy for some fo > and p > and for all s, s S §((i) 

(1.1) vq\s\p < r±(s) : s < + |s|P), 

(1.2) \T^{s)\<v,\l + \sr\ 

(1.3) (^±(,)_^±(^)):(,_^)>0. 

Here we define A : B = tr(AB) for dxd matrices A, B. A typical example is t^(s) = (a^ -|- 

h^\s\'^)~ s with 0=*= > and 6=*= > 0. 

We assume that the velocity field u(x, t) satisfies the following non-Newtonian fluid flow equa- 
tion: 

(1.4) — -Fu-Vn = div(r+(M,n)), divu = on Jl+(t), t > 0, 
at 

dii 

(1.5) — -Fu-Vn = div(r"(M,n)), divu = on f)~(t), t > 0, 

(1.6) n+=n", n • (r+('u,n) - T"(n,n)) = on r(t), t > 0. 

The upper script it in (|1.6p indicates the limiting values approaching to T{t) from f]^(t), respec- 
tively, n is the unit outer normal vector of dVi^{t), H is the mean curvature vector of T{t) and 
> is a constant. The condition ()1.6p represents the force balance with an isotropic surface 
tension effect of the phase boundary. The boundary T{t) is assumed to move with the velocity 
given by 

(1.7) Vt = {u ■ n)n + K2H on T{t), t > 0, 

where /t2 > is a constant. This differs from the conventional kinematic condition (k2 = 0) and is 
motivated by the phase boundary motion with hydrodynamic interaction. The reader is referred 
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to [2T] and the references therein for the relevant physical background. By setting (/? = 1 on Q~^{t), 
(/? = — 1 on Q~{t) and 

r((^,e(n)) = i±^r+(e(^)) + i^r-(e(^)) 
on r2+(i) UQ~{t), the equations (|1.4p - (|1.6|) are expressed in the distributional sense as 

^^g^ ^ + u-Vu = dwT{ip,e{u)) -VU + KiH-W^-^lr^t) onJ]x(0,oo), 

div u = on X (0, oo). 

where T-L'^^^ is the {d — l)-dimensional Hausdorff measure. The expression (jl.Sp makes it evident 
that the phase boundary exerts surface tension force on the fluid wherever H ^ on T(t). Note 
that if T{t) is a boundary of convex domain, the sign of H is taken so that the presence of surface 
tension tends to accelerate the fluid flow inwards in general. We remark that the sufficiently 
smooth solutions of (jl.4p - ()1.7p satisfy the following energy equality, 

(1-9) T-li / \u\'^dx + Kin'^-\T{t))\ = - [ T{ip,e{u)) ■.e{u)dx-KiK2 [ \H\'^ dW^'^ . 
dt [2 Jq } Jn Jr{t) 

This follows from the first variation formula for the surface measure 

(1.10) ^n'^-HT{t)) = - [ Vr -HdW^-^ 

dt Jr(t) 

and by the equations (fLl|) - p^ . 

The aim of the present paper is to prove the time-global existence of the weak solution for 
()1.4p - (jl.7p (see Theorem 12.31 for the precise statement). We construct the approximate solution 
via the Galerkin method and the phase field method. Note that it is not even clear for our 
problem if the phase boundary may stay as a codimension 1 object since a priori irregular flow 
field may tear apart or crumble the phase boundary immediately, with a possibility of developing 
singularities and fine-scale complexities. Even if we set the initial datum to be sufficiently regular, 
the eventual occurrence of singularities of phase boundary or flow field may not be avoided in 
general. To accommodate the presence of singularities of phase boundary, we use the notion of 
varifolds from geometric measure theory. In establishing (|1.7p we adopt the formulation due to 
Brakke [7] where he proved the existence of moving varifolds by mean curvature. We have the 
extra transport effect {u-n)n which is not very regular in the present problem. Typically we would 
only have u G Lfo,([0, oo); W^^P{n)'^). This poses a serious difficulty in modifying Brakke's original 
construction in [7] which is already intricate and involved. Instead we take advantage of the recent 
progress on the understanding on the Allen-Cahn equation with transport term to approximate 
the motion law ()1.7p . 

f +n.V^ = ..(A^-^). (ACT) 

Here W is the equal depth double-well potential and we set W{ip) = (1 — ip'^)'^ /2. When e — )• 0, we 
have proved in [20] that the interface moves according to the velocity p.7p in the sense of Brakke 
with a suitable regularity assumptions on u. To be more precise, we use a regularized version of 
(ACT) as we present later for the result of [20j to be applicable. The result of [20] was built upon 
those of many earlier works, most relevant being [15] which analyzed (ACT) with u = 0, and 

also [laiMiEaisH]. 

Since the literature of two-phase flow is immense and continues to grow rapidly, we mention 
results which are closely related or whose aims point to some time-global existence with general 
initial data. In the case without surface tension (ki = H2 = 0), Solonnikov [31] proved the 
time-local existence of classical solution. The time-local existence of weak solution was proved by 
Solonnikov [32], Beale [S|, Abels P], and others. For time-global existence of weak solution, Beale 
[6] proved in the case that the initial data is small. Nouri-Poupaud [26] considered the case of 
multi-phase fluid. Giga-Takahashi [11] considered the problem within the framework of level set 
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method. When ki > 0, K2 = 0, Plotnikov [27] proved the time-global existence of varifold solution 
for d = 2, p > 2, and Abels [2J proved the time-global existence of measure-valued solution for 
d = 2,3, p > ■j^. When Ki > 0, K2 > 0, Maekawa [22j proved the time- local existence of classical 
solution with p = 2 (Navier-Stokes and Stokes) and for all dimension. Abels-Roger [3] considered a 
coupled problem of Navier-Stokes and Mullins-Sekerka (instead of motion by mean curvature in the 
present paper) and proved the existence of weak solutions. As for related phase field approximations 
of sharp interface model which we adopt in this paper, Liu and Walkington [21] considered the 
case of fluids containing visco-hyperelastic particles. Perhaps the most closely related work to the 
present paper is that of Mugnai and Roger [25| which studied the identical problem with p = 2 
(linear viscosity case) and d = 2,3. There they introduced the notion of velocity and showed 
that ()1.7p is satisfied in a weak sense different from that of Brakke for the limiting interface. Kim- 
Consiglieri-Rodrigues [16] dealt with a coupling of Cahn-Hilliard and Navier-Stokes equations to 
describe the flow of non-Newtonian two-phase fluid with phase transitions. Soner [33] dealt with 
a coupling of Allen-Cahn and heat equations to approximate the Mullins-Sekerka problem with 
kinetic undercooling. Soner's work is closely related in that he showed the surface energy density 
bound which is also essential in the present problem. 

The organization of this paper is as follows. In Section 2, we summarize the basic notations and 
main results. In Section 3 we construct a sequence of approximating solutions for the two-phase 
flow problem. Section 4 describes the result of PO] which establishes the upper density ratio bound 
for surface energy and which proves (|1.7p . In the last Section 5 we combine the results from Section 
3 and 4 and obtain the desired weak solution for the two-phase flow problem. 

2. Preliminaries and Main results 

For dx d matrices A, B we denote A : B = ti (AB) and \ A\ := V A : A. For a G M'^, we denote 
by a (8) a the d x d matrix with the i-th row and j-th column entry equal to OjOj. 

2.1. Function spaces. Set Q, = T'^ throughout this paper. We set function spaces for p > 
as follows: 

V = jz; G C^in^; dwv = o} , 

for s G Z+ U {0}, W'^iQ) = {v : V^v G LP{n) for < j < s}, 
V''P = closure of V in the H^"'P(f])'^-norm. 

We denote the dual space of V^'^ by (^'P)*. The inner product is denoted by (•,•). Let xa be 
the characteristic function of A, and let | Vxa| be the total variation measure of the distributional 
derivative Vxa- 

2.2. Varifold notations. We recall some notions from geometric measure theory and refer to 
[U [71 [30] for more details. A general k-varifold in R"' is a Radon measure on R*^ x G{d, k), where 
G{d, k) is the space of A;-dimensional subspaces in R''. We denote the set of all general /c-varifolds 
by Vfc(R'^). When 5 is a A;-dimensional subspace, we also use S to denote the orthogonal projection 
matrix corresponding to R'^ — )• S. The first variation of V can be written as 

SV{g)= [ Vgix):SdVix,S) = - [ g{x) ■ H{x) d\\V\\ix) if || « ||y||. 

Here V G Vfc(R'^), \\V\\ is the mass measure of y, g G C^{W^)'^, H = Hy is the generalized mean 
curvature vector if it exists and \\5V\\ <^ \\V\\ denotes that \\6V\\ is absolutely continuous with 
respect to \\V\\. 

We call a Radon measure /i k-integral if // is represented as /i = 8'H''[x, where X is a countably 
fc-rectifiable, T^'^-measurable set, and 6 G Ll^^{'H^[x) is positive and integer-valued Ti^ a.e on X. 
T-L^^x denotes the restriction of T-L^ to the set X. We denote the set of fc-integral Radon measures 
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by TM.}^. We say that a /c-integral varifold is of unit density if ^ = 1 l-O' a.e. on X. For each such 
fc-integral measure ^ corresponds a unique fc-varifold V defined by 



ct){x, S) dV{x, S) = / T^^u) dfi{x) for </> € Cc(K'' x A;)), 

'xG{d,k) JK.'i 

where T^fj, is the approximate tangent fc-plane. Note that /i = ||1^||. We make such identification 
in the following. For this reason we define Hf^ as Hy (or simply H) if the latter exists. When X 
is a submanifold without boundary and 6 is constant on X, H corresponds to the usual mean 
curvature vector for X. In the following we suitably adopt the above notions on Q = T'^ such as 
Vfc(r2), which present no essential difficulties. 

2.3. Weak formulation of free boundary motion. For sufficiently smooth surface r(t) moving 
by the velocity ()1.7p . the following holds for any G C'^{Q;'R'^) due to the first variation formula 
(fTTnl) : 



(2.1) ^/ (^dn'^-^ < [ {-(l)H + V^)-{K2H+{u-n)n}dn'' 

dt Jr(t) Jr(t) 



id-l 

'm J^it) ' 

One can check that having this inequality for any (p G C^(r2;M^) implies (|1.7p thus ()2.ip is 
equivalent to ()1.7p . Such use of non-negative test functions to characterize the motion law is due 
to Brakke [7j where he developed the theory of varifolds moving by the mean curvature. Here we 
suitably modify Brakke's approach to incorporate the transport term u. To do this we recall 

Theorem 2.1. ( Meyers- Ziemer inequality) For any Radon measure fi on W^with 

D = sup j—^ < oo, 

r>0,a;GR'* ^d-lf 

we have 



(2.2) / \(j)\dfi <cD \V(t)\dx 



for (j) G C^(M ). Here c depends only on d. 

See j24] and |35|, p. 266]. By localizing (|2.2p to = T*^ we obtain (with r in the definition of D 
above replaced by < r < 1/2) 



(2.3) ^ \ct>? di, < cD [ml^^n) + ||V0||i.(^)) 



where the constant c may be different due to the localization but depends only on d. The inequality 
()2.3p allows us to define d/j, for (p G Ty^'^(ri) by the standard density argument when D < oo. 

We define for any Radon measure /x, u G L^(O)'^ and (/> G C^(0 : M^) 

(2.4) B{n, u, (/))= / {-(t>H + V0) • {k2H + (n • n)n} d^ 

Jn 

if /i G with generalized mean curvature H G L'^{fi)'^ and with 

(2.5) sup -T—r < oo 

i>r>0,xeC^d-l^ 

and u G W^'"^ {Q)'^ . Due to the definition of IAid^i{^), the unit normal vector n is uniquely 
defined n a.e. on 0, modulo it sign. Since we have {u,n)n in (|2.4p . the choice of sign does not 
affect the definition. The right-hand side of (|2.4p gives a well-defined finite value due to the stated 
conditions and (j2.3p . If any one of the conditions is not satisfied, we define B{fi, u, (p) = — oo. 
Next we note 
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Proposition 2.2. For any < T < oo and p > 

|nGL?^([0,r];yi'f); ^ GL^([0,r];(yi'f)*)| -^C([0,T]; ^0,2). 

The Sobolev embedding gives V^'"^ ^ V^''^ for such p and we may apply [23^ p. 35, Lemma 
2.45] to obtain the above embedding. Indeed, we only need p > for Proposition 12.21 to be 
and we have > -j^- Thus for this class of u we may define u{-,t) G V^'^ for all t G [0,T] 
instead of a.e. t and we may tacitly assume that we redefine u in this way for all t. For {fJ't}te[o,oo)^ 

u G Lf^^{[0,oo);V^'P) with ^ G ([0, oo); (F^'f)*) for p > and (p G C^{n;R+), we define 
B{fit, u{-,t), (p) as in (glD for all t > 0. 

2.4. The main results. Our main results are the following. 

Theorem 2.3. Let d = 2 or 3 and p > ^ . Let ft = T'^. A ssume that locally Lipschitz functions 

: S{d) -)> S(d) satisfy (frT]l - (fL3l) . For any initial data uq G and n+{0) C /iamng 
boundary d^l~^{0), there exist 

(a) 7xGL-([0,oo);yO-2)nLP^([0,oo);yi'?') ^zi/i ^ g ([0, oo); (y^'f)*), 

(b) a family of Radon measures {fJ't}te[o,oo) '^^^^ £ /or a.e. t G [0,oo) and 

(c) V9 G X [0, oo)) n L°°([0, oo); SF(l^)) n C;t([0, oo); Li(f])) 
suc/i i/iat t/ie following properties hold: 

(i) T/ze triplet {u{-,t), ^p{-,t), fJ't)t£[o,oD) weak solution of (|1.8|) . More precisely, for any 
T > we have 

(2.6) I f —u ■ ^ + {u -Vu) • V + T{ip,e{u)) : e{v) dxdt = f UQ-v{0)dx+ f f KiH-vdfj-tdt 

Jo Jn Jn Jo Jn 

for any v G C~([0,T];V) such that v{T) = 0. Here H G L^ {[0, oo); L'^iptf) is the 
generalized mean curvature vector corresponding to p,f 

(ii) The triplet {u{-,t), ip{-,t), /Ut)te[o,oo) satisfies the energy inequality 

- I \u{-,T)\'^ dx + KifiTi^) + I I T{ip,e{u)) : e{u) dxdt + KiK2 [ [ \H\'^ dfifdt 
2 Jn Jo Jn Jo Jn 



(2.7) 



<l [ luol"^ dx + nin'^-\dn+ (0)) =:Eq 
^ Jq 



for all T < oo. 
(in) For all < ti < t2 < oo and (j) G C^(ri;R+) we have 

(2.8) ^it, (0) - lit, ((/>)< / ' S(/Xi, u{-,t), (/>) dt. 



ti 



Moreover, B{fit, u{-,t), 0) G ^^^^([0, oo)). 

(iv) We set Dq = supq<^<]^/2,3:go — ^^1-^'^^^''^ ■ For any < T < oo, there exists a 
constant D = D{Eq, Do,T,p, uq, ki, K2) such that 

sup ^^^^^ < D 

0<r<l/2, xGQ ^d-lf 

for all t G [0,r]. 

(v) The function ip satisfies the following properties. 

(1) = ±1 a.e. on for all t G [0, 00) . 

(2) <p{x,Q) = Xq.+ {0) - Xn\n+{0) a.e. on n. 

(3) spt|Vx{<^(.,t)=i}| C spt//t for all t G [0,oo). 
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(vi) There exists 

Ti = Ti{Eo,Do,p,i^o,Ki,K2) > 
such that /it is of unit density for a.e. t G [0, Ti]. In addition |Vx{,^(.,t)=i}| = /it for a.e. 

te[o,Ti]. 

Remark 2.4. Somewhat different from u = case we do not expect that 
(2.9) hmsup — <B{iit,u{-,t),(j)) 

holds for all t > and 4> £ C^(il;IR+) in general. While we know that the right-hand side is < oo 
(by definition) for all t, we do not know in general if the left-hand side is < oo. One may even 
expect that at a time when \'Vu{-,t)\^ dx = oo, it may be oo. Thus we may need to define (|1.7p 
in the integral form (|2.8p for the definition of Brakke 's flow. Note that in case u = 0, one can 
show that the left-hand side of (|2.9p is < oo for all t >0 (see [7\). 

Remark 2.5. The difficulty of multiplicities have been often encountered in the measure-theoretic 
setting like ours. Varifold solutions constructed by Brakke [7J have the same properties in this 
regard. On the other hand, (vi) says that there is no 'folding' for some initial time interval [0,Ti] 
at least. 

Remark 2.6. In the following we set ki = ^2 = 1 for notational simplicity, while all the argument 
can be modified with any positive ki and K2 with no essential differences. On the other hand, 
their being positive plays an essential role, and most of the estimates and claims deteriorate as 
Ki, — 7- and fail in the limit. How severely they fail in the limit may be of independent interest 
which we do not pursue in the present paper. Note that k,2 = limit should correspond precisely 
to the setting of Plotnikov [27\ for d = 2. 

We use the following theorem. See [23, p. 196] and the reference therein. 

Theorem 2.7. (Korn's inequality) Let 1 < p < oo. Then there exists a constant ck = c{p,d) 
such that 

\\v\\w^,P^n) ^ Cit(||e(f)||^^(^) + ||^;||^i(f^)) 

holds for all v G W^^P{n f. 

3. Existence of approximate solution 

In this section we construct a sequence of approximate solutions of (jl.4p - (|1.7p by the Galerkin 
method and the phase field method. The proof is a suitable modification of [18j for the non- 
Newtonian setting even though we need to incorporate a suitable smoothing of the interaction 
terms. 

First we prepare a few definitions. We fix a sequence {sj}^]^ with limj^oo £i = and fix a 
radially symmetric non- negative function G C^{W^) with spt^ C -Bi(O) and f ^ dx = 1. For a 
fixed < 7 < ^ we define 

(3.1) 

We defined C» so that f dx = 1, < c{d)e~^ and |VCH < c{d)e';^~^''^ . 

For a given initial data fi+(0) C VL with boundary 0^7+ (0), we can approximate ri+(0) in C 
topology by a sequence of domains r2*"'~(0) with boundaries. Let d^{x) be the signed distance 
function to 90*+ (0) so that d\x) > on O*+(0) and d*(x) < on il*-(0). Choose 5* > so that 

is function on the ^^-neighborhood of 9$7*"'"(0). Now we associate {ej}^^ with r2*"'"(0) by re- 
labeling the index if necessary so that limj_^oo Ei/^' = and limj_>.oo ei~^\V^d^\ = for j = 2, 3 on 
the 6*-neighborhood of dQ'^~^{0). Let h G C°°(M) be a function such that h is monotone increasing, 
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h{s) = s for < s < 1/4 and h{s) = 1/2 for 1/2 < s, and define h{-s) = -h{s) for s < 0. Then 
define 

(3.2) if^Q^ix) = ta.nh{b' h{d'{x)/b')/ei). 

Note that we have (/?q* G C^(0) and el\V^ipy-\ for j = 1, 2, 3 are bounded uniformly independent 
of i. The well-known property of phase field approximation shows that 

(3.3) hm - ixnno) " xn-{o))llLi{Q) = 0, - (f^^^ + dx ^ n'-'ien+m 

as Radon measures. Here a = J_^' ./2W{s)ds. 

For V'^''^ with s > | + 1 let {o;*}^^ be a set of basis for V^''^ such that it is orthonormal in V^'"^. 
The choice of s is made so that the Sobolev embedding theorem implies W^~^''^{Q) L°°(J7) thus 

Let Pj : V^'"^ — )• V^'^'^ = span{c<Ji, W2, • • • ,0Ji} be the orthogonal projection. We then project the 
problem ()1.4p - p.7p to ' by utilizing the orthogonality in V^''^. Note that just as in [T8], we 
approximate the mean curvature term in (|1.8p by the appropriate phase field approximation. We 
consider the following problem: 

(3.4) -^ = Pi [dwT{ip^\e{u^')) -u""' ■ Vu^' - -^div {{Vip^' Vip""') * on Q x [0,oo), 

(3.5) u''{-,t) eV^'^ fort>0, 



(3.6) + {u^^ Vv9^' = Aif^^ onQx [0, oo), 



ot e.^ 
(3.7) w^'(x,0) = Pi'uo(x), (/?^'(x,0) = v3o'(x) on 17 

Here * is the usual convolution. We prove the following theorem. 

Theorem 3.1. For any i € N, uq G V^'"^ and ip'^ , there exists a weak solution {u''%p^') of ([3^ 
(15:71) such thatu'^ G L°°([0,oo);F0'2) nLf^^([0,oo);yi'P), < 1, p"' G L°°([0, oo); C3(J7)) and 

¥GL-([0,oo);C7Hf^)). 

We write the above system in terms of ti^' = Yl\=i (i)wfe(x) first. Since 

i 

(u^« • Vn=% Wj) = 4'(t)cf (t)(wfc • Vw,, ojj), 



kl=l 



ei(div ((Vv?"' (g) V^J^O * ^j) = -^i / (Vv?"' ® Vy?^') * C ■ Vcjj drc, 

(divr((^^% e(M^')), Wj) = - / T{p'^\e{u'^^)) : e{ujj) dx 

Jn 

for j = 1, • • • , i, ()3.4p is equivalent to 



dt ^ 

(3.8) 



d f 

cf{t) = - T{p'^,e{u'^)) : e{ujj)dx - ^ cl'{t^'{t){ujk ■ Vui, ujj] 

-'^ kl=l 



+ ^ / (Vv?^' Vp^^) * : dx = A'Mt) + BkijclUt^^ (t) + Df{t). 
Jn 



Moreover, the initial condition of c^* is 



c^'(O) = (no, Wj) for j = 1,2, 
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We also set 



Eo = n'^-\dn+io)) + l I \uo\^dx 

^ Jn 



and note that 
(3.9) 

by p.3p and by the projection Pi being orthonormal. 
We use the following lemma to prove Theorem 13.11 

Lemma 3.2. There exists a constant Tq = ro(£'o, i, t'cp) > such that ()3.4p - ()3.7p with ()3.9p has 
aweak solution mOx[0,To] such thatu^^ G L~([0, Tq]; F°'2)nLP([0, Tq]; y^'f), < 1, 

if'' G L~([0,ro];C3(0)) and ^ E L^i[0,To];C\n)). 

Proof. Assume that we are given a function u{x,t) = X]j=i c^'(t)cjj(x) G C^^'^{[0,T];V^''^) with 

/ . \ 1/2 . 

(3.10) c^^»(0) = (no, coj), max i fl 1^'^*)!' + E n/^^' < 
We let (/9(x,t) be the solution of the following parabolic equation: 

(3.11) dt 



^{x,0) = ifl'ix). 

The existence of such ip with ly?! < 1 is guaranteed by the standard theory of parabolic equations 
([17j). By (j3.1ip and the Cauchy-Schwarz inequality, we can estimate 



Since for any t € [0, T] 

i 

ll^*C^1li-(n) ^ ^ ^^^r^^^ max ||tjj(x)|||^(f^)^|c^*(t)P < c{i)Eo, 



i=i 



This gives 



0<t<T f V 2 / 



Hence as long as T < 1, 



(3.13) \Dnt)\ < cWVoojUoo^n)- [ f ei\V^{y)fC'{x - y) dydx < c(i)e'^«^»Eo 



^ Jn Jn 

by Vcoj e L°°(fi)'^' and (l3T2]l . 

Next we substitute the above solution 99 into the place of ip^% and solve p.Sp with the initial 
condition c^'(O) = (uq, ujj). Since r is locally Lipschitz with respect to e(u), there is at least 
some short time Ti such that (j3.8p has a unique solution c^'(i) on [0,Ti] with the initial condition 
c^'(O) = (uo) '^i) for 1 < j < i. We show that the solution exists up to Tq = To(i, E'ojPj ^0) 
satisfying (IXTUT) . Let c(t) = \ Y!j=i Then, 
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By ([LID ^"^c^j < hence 

±c{t)<c{^,E,){?^^/^+c^'\ 

Therefore, 

(3.14) arctan c{t) < arctan \/Eq + 2c(i, EQ)t. 

We can also estimate \dc^j^ /dt\ due to ()3.8p . (jS.lSp . ()3.14p and ()1.2p depending only on EQ,i,p,UQ. 
Thus, by choosing Tq small depending only on EQ,i,p,VQ we have the existence of solution for 
t G [0,To] satisfying ()3.10p . We then prove the existence of a weak solution on x [0,To] by using 
Leray-Schauder fixed point theorem (see [E]). We define 

i 

u{x,t) = ^c^'(t)cjj(x) 

J=l 

and we define a map £ : u i— )• tt as in the above procedure. Let 

1/2 



V{To) := l u{x,t) = ^Cj{t)ujj{x) ; max \l^\cj{t)('' 



0<4i<t2<Top^ 1^1 -^2^'^ ' 



Then V{Tq) is a closed, convex subset of Ci/2Qo,ro]; V;°'^) equipped with the norm 

||n||v.(To)= max -y|c,(t)n + sup y hpl^^J^^ 

and by the above argument C : V(Tq) — )• V{Tq). Moreover by the Ascoli-Arzela compactness 
theorem £ is a compact operator. Therefore by using the Leray-Schauder fixed point theoremC£ 
has a fixed point n^* E ^(^o)- We denote by cp'^^ the solution of (|3.6p and (|3.7p . Then {u'^\if 
is a weak solution of (j3.4p - (|3.7p in ^2 x [0,ro]. Note that we have the required regularities for if 
due to the regularity of * C,^^ in x and by the standard parabolic regularity theory. □ 

Theorem 3.3. Let (u^^cp^') be the weak solution of ([33]) -([321 with <^M) in Q x [0,T]. Then the 
following energy estimate holds: 



(3.15) 



+ ^ ^^(^Ac^^'- ^y^ ^ +i/o|e(u^Or^^^'^* <^o + o(l). 



Moreover for any < Ti < T2 < oo 



(3.16) / ' h''(-,t)Ci,P(n) < Ci^{^o"'^o + (T2 - ri)4} + 0(1) 
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Proof. Since {u^',ip^^) is the weak solution of (|3.4p - (j3.7p . we derive 



(3.17) 



dt Jo cr \ 2 £i J 2 



dx 



+ j |divr(99^%e(n''')) - u"' ■ Vu"' - — div((V(/?^' (g) Vip"') * C')^ ■ u"' dx = h + h- 
Since div (n^* * C^*) = (divu"^') * = 0, 

ah = - [ Ei (aip"' - I^IM^ dx + Si / (n^' * C') • V(^^'A(^^' dx. 
For /2, with (fmi 

divT(99'^%e(u'^^)) • u^' dx = - / t((^^\ e(ti'^')) : e(ti'^') f^^; < -z^o / |e(ti'^')r t^^;. 
Moreover the second term of I2 vanishes by divu"^' = and 

- J ejdiv {Vif^' (g) V(p^' ■ -u"' dx = - J Si ^V ^^'^ ^ + V(/3^* A(/7^"^ * • u"' dx 

= -Ei / (u^" * CO • V(/5^'A(/9^' dx. 
Hence (|3.17|) becomes 

' ( *i ^ I + ) + dx<- I -^{Aip'^ +z^o|e(u^Ol^c^a;. 



dt Jn(J \ 2 Ei J 2 V , 

Integrating with respect to t over t G [0,T] and by (|3.9|) . we obtain (|3.15p . The proof of (|3.16p 
follows from (|3.15p and Theorem 12. 7[ □ 

Proof of Theorem \3.1\. For each fixed i we have a short time existence for [0, Tq] where Tq 
depends only on i,EQ,p,h'Q at t = 0. By Lemma [3.31 the energy at t = Tq is again bounded by 
Eq + 0(1). By repeatedly using Lemma 13.21 Theorem 13.11 follows. □ 



4. Proof of main theorem 

In this section we first prove that {99^'}^]^ in Section 3 and the associated surface energy 
measures {fil'}°2^^ converge subsequentially to if and fit which satisfy the properties described in 
Theorem 12.31 Most of the technical and essential ingredients have been proved in [20j and we 
only need to check the conditions to apply the results. We then prove that the limit velocity field 
satisfies the weak non-Newtonian flow equation, concluding the proof of Theorem 12.31 
First we recall the upper density ratio bound of the surface energy. 

Theorem 4.1. Theorem 3.1]^ Suppose d>2,n = T'^,p> i > 7 > 0, 1 > e > and 

if satisfies 

(4.1) ^ + u-Vip = Aif-^-^ onl7x[0,r], 

(4.2) (/9(x,0) = ^o{x) on 0, 
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where V'u, V^(p,V^(pt G C{n x [0,r]) for < i, k < I and < j < 3. Let fit be the Radon 
measure on $7 defined by 

,4.3) = + 



n cr V 2 e 

for (f) £ C{n), where a = J\ ^2W{s)ds. We assume also that 
(4.4) sup \ ipo\ < 1 and supe*|VVo| < ci for 1 < i < 3, 

,4.5) " ,,pf^_M<,-., 



n 



2 e 



(4.6) sup {e^|u|, e^+^|Vn|} < C2, 

Ox[0,T] 

(4.7) /" ||n(-,t)||^;^,,,,(^)(it <C3. 







Define for t G [0, T] 



(4.8) D(t)=maxJ sup ^/it(5,(x)), 1 I , I)(0) < Dq- 

T/ien i/iere exisi ei > which depends only on d, p, W, c\, C2, C3, Dq, 7 and T, and C4 which 
depends only on C3, d, p, Dq and T such that for all < e < ei, 

(4.9) sup D{t) < C4. 

0<t<T 

Using this we prove 

Proposition 4.2. For {(^^'}^| in Theorem \ 3.1l define as in (|4.3p replacing ip by ip^'- , and 
define D^'{t) as in ()4.8p replacing fit by fif'. Given <T < 00, there exists C5 which depends only 
on Eq, vq, 7, ^5 d, p and W such that 

(4.10) sup D^^(i) < C5 

0<t<T 

/or a// sufficiently large i. 

Proof. We only need to check the conditions of Theorem 14.11 for (p'^^ and Note that u in 
gH) is replaced by u^' * C*"'- We have d>2,n = T'^,p> ^ > 7 > 0, 1 > e > and (gl]) and 
()4.2p . The regularity of functions is guaranteed in Theorem 13.11 With an appropriate choice of ci , 
()4.4p is satisfied for all sufficiently large i due to the choice of in ()3.2p . The sup bound ()4.5p is 
satisfied with even on the right-hand side instead of £~"' . The bound for u^* * C*"* (|4.6p is satisfied 
due to (|3.ip and (|3.15p . and ()4.7p is satisfied due to ()3.16p . Thus we have all the conditions, and 
Theorem 14.11 proves the claim. □ 

We next prove 

Proposition 4.3. For {u^^ * C^^}iZi ^'^ Theorem \ 3.1l there exist a subsequence (denoted by the 
same index) and the limit u G L°°([0, oo)\ V^'"^) n L^^^([0, 00); y^'^) such that for any < T < 00 
(4.11) 

ys, ^Qs, weakly in L^{[Q,T\;W^^'P{9.)'^), u"' * u strongly in L^{[0,T]; L^{nf). 

Proof. Let ip G V'"^ with ||V'||v''.2 < 1- With ([331), ([33]) and integration by parts, we have 

T((y9^\e(ti^')) + —{Vif^' V(/?^') * C',VPitp) . 

a / 
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Here we remark that 

||VPiVllL°°(f7) < ciMPii^Ww'^n) < c(d)||V'||iy=>2(Q) = c(d)||^/;||y«,2 < c{d) 
by s > ^ and properties of Pi (see [TO] or [Ml p.290]). Thus by (USD and (fXT^ . we obtain 

' ^u'^'■ 



,tp] < c{d,p,uo) (I + E0 + \\u 



Again using (|3.15|) and integrating in time we obtain 

p 

p-i 



(4.12) 



f 




Jo 


dt 



dt < c{d,p,Eo,iyo,T). 



Now we use Aubin-Lions compactness Theorem [191 p. 57] with Bq = V^'"^, B = V'^''^ C L?{Q)'^, 
Bi = (y*'^)*, po = p and pi = Then there exists a subsequence still denoted by {u^^}'^^ such 
that 

u'^^u in LP{[0,T];L'^{ny). 

Since we have uniform L°°{[0,T]; L'^{Q)'^) bound for the strong convergence also holds in 
L'^{[0,T];L^{ny). Note that we also have proper norm bounds to extract weakly convergent 
subsequences due to ()3.15p . For each r„ which diverges to 00 as n — )• 00, we choose a subsequence 
and by choosing a diagonal subsequence, we obtain the convergent subsequence with ()4.1ip with 
ti^' instead of u^* * C^' • It is not difficult to show at this point that the same convergence results 
hold for u^' * C as in (fiTT]l . □ 

Proof of main theorem. At this point, the rest of the proof concerning the existence of the 
limit Radon measure and the limit if = limj_j.oo V'^' and their respective properties described in 
Theorem 12 . 3 1 can be proved by almost line by line identical argument in |20[ Section 4, 5]. The only 
difference is that the energy Eq in [20j depends also on T, while in this paper Eq depends only on the 
initial data due to (|3.15|) . This allows us to have time-global estimates such as u E L°°{[0, 00); V^'^) 
and if G L°°([0, 00); BV{i})). The argument in [2Uj then complete the existence proof of Theorem 
12.31 (b). (c) along with (iii)-(vi). We still need to prove (a), (i) and (ii). 

Due to (j4.12p . p.2p and (j3.16p we may extract a further subsequence so that 
(4.13) 

^''''^^^ L—{[{),T]- {V'^Y), T{ip'\e{u'^)) ^ f weakly in L—{[{),T]-L—{^f). 
For ujj £ {j = !,■■■) and hG C^{{0,T)) we have 

[ diY{{\7(p^' (g)Vip''')*C')-hujjdx= [ ( Aip"' _ ^ ^ \ Vy?^' • hujj * dx 
Jn Jn\ J 

by integration by parts and diw uoj = 0. Thus the argument in [19', p. 212] and the similar conver- 
gence argument in |20j show 



(4.14) j J {u -Vu) ■ huij + hf : e{ujj)dx^ dt = J J H ■ htoj dfitdt. 



Again by the similar argument using the density ratio bound and Theorem 12.11 one show by the 

du 

at 



density argument and that ^ G Lp-^{[0,T]; (F^'P)*) and 



(4.15) j |(^^'^^+ j {u-'^u)-v + f ■.e{v)dx^dt = j j H-vdfitdt. 
for ah V G LP{[0,T];V^'P). We next prove 

(4.16) [ [ T:e{v)dxdt= [ [ T{if,e{u)) : e{v) dxdt 

Jo Jn Jo Jn 
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for all V e C^{{0,T);V). As in [H p.213 (5.43)], we may deduce that 

(4.17) \\Hh)\\l2^n) I' IJ'- - r In^ -udixtdt + \\Hml2in) 
for a.e. ti G [0,T]. We set for any v G V^'^ 

(4.18) r / (T(<^^%e(^z^O)-r(¥^'%e(t'))):(e(n-0-e(z;))(ixdt + i||^x^>(ti)||i.(^). 

"/ J 

The property (jl.Sp of e(-) shows that the first term of ()4.18p is non-negative. We may further 
assume that it'^'(ii) converges weakly to u{ti) in L^(ri)'^ thus we have 

(4.19) linimf4i>i||tx(ti)||i,(^). 
By (|3.4p we have 



111 /„N ii9 £i 



4^ =2 ll^'"(0)lli^(f^) - - / / div((V(^^» ^ Vv^^O * CO • 

/ r(99^% e(u^^)) : e(t;) + r(99'^% e(7;)) : {e{v!^^) - e{v)) dxdt 
./f7 

which converges to 

(4.20) A^^ = \\\u{0)\\l2,^.+ r [ H-udntdt- [' I f : e{v) + T{^,e{v)) : {e{u) - e{v)) dxdt. 
2 ^ ' . n ./o ./n ./o 



'0 Jn Jo Jn 

Here we used that f'^' converges to 99 a.e. on 17 x [0,T]. By (|4.17p . (|4.19p and (|4.20p . we deduce 
that 

/ / (t - T{^p, e{v))) : (e(n) - e{v)) dxdt > 0. 
Jo Jn 

By choosing v = u + ev, divide by e and letting e — )• 0, we prove ()4.16p . Finally, ()2.7p follows from 
()4.16p . the strong L"^(Ox [0, T]) convergence of ip''\ the lower semicontinuity of the mean curvature 
square term (see [20]) and the energy equality appearing in Theorem 13.31 This concludes the proof 
of Theorem O □ 
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